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Abstract 

Downlink scheduling schemes are well-known and widely investigated under the assumption that 
the channel state is perfectly known to the scheduler In the multiuser MIMO (broadcast) case, downlink 
scheduling in the presence of non-perfect channel state information (CSI) is only scantly treated. In this 
paper we provide a general framework that addresses the problem systematically. Also, we illuminate the 
key role played by the channel state prediction error: our scheme treats in a fundamentally different way 
users with small channel prediction error ("predictable" users) and users with large channel prediction 
error ("non-predictable" users), and can be interpreted as a near-optimal opportunistic time-sharing 
strategy between MIMO downlink beamforming to predictable users and space-time coding to non- 
predictable users. Our results, based on a realistic MIMO channel model used in 3GPP standardization, 
show that the proposed algorithms can significantly outperform a conventional "mismatched" scheduling 
scheme that treats the available CSI as if it was perfect. 
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I. Introduction 

Under perfect knowledge of the downlink channels, the resource allocation problem in a 
Multiuser MIMO (MU-MIMO) downlink has been widely investigated under various precoding 
and beamfomring schemes (see for example [1], [2], [3], [4], [5], [6], [7], and references therein). 
In practice, the Channel State Information (CSI) is obtained through some form of training 
and feedback. In Time-Division Duplexing (TDD) systems, the base station (BS) can learn the 
downlink channel coefficients in "open-loop" mode, by exploiting the uplink pilot symbols and 
channel reciprocity (e.g., [8], [9]). In Frequency-Division Duplexing (FDD) systems, since uplink 
and downlink take place in widely separated frequency bands, the downlink channel coefficients 
must be learned in "closed loop" mode, via some explicit CSI feedback scheme (e.g., [10], [11] 
and references therein). In both cases, the CSI available to the BS can be seen as some sort of 
"noisy" version of the true channel coefficients. 

The key impact of CSI quality on the performance of MU-MIMO is evidenced in the relevant 
regime of medium-to-high SNR. In [10] it is shown that the gap between the sum capacity under 
perfect CSI and the sum-rate achievable by a simple linear Zero-Forcing Beamforming (ZFBF) 
scheme with non-perfect CSI takes on the form AR = min{M, log(l + KiCTgSNR) + K2, 
where ki, K2 are constants that depend on the particular CSI training and feedback scheme used, 
and where erf denotes the Mean-Square Error (MSE) between the true channel coefficients and 
the CSI available to the BS. Since the sum-rate of a MU-MIMO downlink channel in high-SNR 
behaves like min{M, K} log(l + kSNR), for some constant k, it follows that must decrease 
at least as fast as SNR~^ in order to preserve the optimal 0(min{M, .fr} logSNR) increase of 
the sum-rate with SNR. Under both TDD and FDD, the main source of CSI estimation error 
consists of the delay introduced by the estimation/feedback scheme in the presence of time- 
varying wireless channels [8], [10], [11], [12], [13]: even after neglecting all other sources of 
suboptimality, such as channel state quantization, feedback errors, and so on, the MSE cannot 
be less than the channel prediction error from noisy pilot symbols. This estimation-theoretic 
quantity represents a fundamental lower bound to the accuracy of CSI. 

We performed an extensive study of MIMO channel prediction based on the 3GPP Spatial 
Channel Model given in [14]. Our results, summarized in Section [V-Bi show that channel pre- 
diction is generally quite accurate with the exception of a specific class of channels characterized 
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by a large Doppler spread (high user mobility) and clustered angles of arrival. For such channels, 
the channel prediction MSB is very large, no matter which prediction method is used, as reflected 
by a Cramer-Rao bound analysis not included in this work for the sake of space limitation [15]. 
These results suggest that users can be classified according to their channel prediction MSB, and 
that this classification effectively reduces to only two extreme classes of "non-predictable" (high- 
mobility and clustered angle of arrivals) and "predictable" (all other cases). This simplification is 
instrumental to the main contribution of this paper: a simple and efficient MU-MIMO downlink 
scheduling scheme that takes explicitly into account the CSI quality. 

Downlink scheduling aims at making the system operate at a desired point on the ergodic 
(or long-term average) achievable rate region of the system, for a given physical layer signaling 
scheme. The operating point reflects some form of "fairness," corresponding the maximization 
of a concave non-decreasing utility function of the ergodic rates. Although a direct maximization 
is typically hopelessly complicated, the optimal point is implicitly achieved using a stochastic 
optimization approach [16], [17], [18]. We solve this problem for the case of MU-MIMO with 
non-perfect CSI in Section [nil Then, based on the general solution, we find a practical simplified 
scheduling policy under the assumption, motivated before, that the users can be partitioned into 
two classes with either small or large channel prediction MSB. The resulting scheduling algorithm 
can be regarded as an opportunistic MIMO "multi-mode" scheme that selects at each scheduling 
slot either a MU-MIMO downlink beamforming mode that performs spatial multiplexing to a 
subset of predictable users, or a single-user space-time coding mode that serves a single selected 
non-predictable user. 

With respect to existing literature, we notice that downlink scheduling with non-perfect CSI 
has been treated mainly in the case where all users have the same CSI quality. Static mode- 
switching criteria have been studied for example in [19], [13] where the number of users to be 
simultaneously served is optimized depending on the CSI quality and channel SNR. In contrast, 
the present work presents a dynamic scheduling policy that can handle users with very different 
CSI qualities at the same time, and allocates opportunistically the signaling modes (namely: 
spatial multiplexing and space-time coding) over time and across the users. The fundamental 
role of channel prediction in downlink scheduling schemes was noticed before, e.g., in [13], 
[12]. In particular, [12] proposes a channel-predictive proportional fair scheduling rule, without 
analytical proof, for the scalar (not MIMO) case. In comparison with these works, here we 
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provide a general framework for downlink scheduling with non-perfect CSI that applies to MU- 
MIMO and to a wide class of fairness utility functions. Also, we present novel rigorous results 
on system stability and performance bounds of the proposed scheduling algorithms. 

Numerical results are provided for two relevant fairness utility functions reflecting proportional 
fairness and max-min fairness (referred to as "hard-fairness"). It should be noticed, though, that 
our framework can be applied to any concave non-decreasing utility function. Results based 
on a realistic channel model [14] and actual channel state prediction algorithms (see details in 
Section jV]) show that the proposed approach achieves very significant improvement with respect 
to a conventional mismatched scheme that treats the available CSI as if it was perfect. 

II. System set-up 

We consider a MIMO downlink channel with a BS equipped with M antennas and K single- 
antenna UTs. The channel is assumed frequency flalQ and constant over "slots" of length T ^ 1 
symbols (block-fading model). The received complex baseband discrete-time signal at the fc-th 
UT during block t is described by 

|/fc,,(t) =h^(t)Xi(t) + 2fc,i(t), 2 = 1,...,T (1) 

where ^ denotes Hermitian transpose, t tics at the slot rate, i tics at the symbol rate, k denotes 
the user index, hfc(t) G C^''^ is the channel vector from the BS antenna array to the A;-th receiver 
antenna, Xj(t) G C^^ is the transmit signal vector transmitted at symbol interval i of slot t, and 
Zk,i{t) ~ 6^(0, A'^o) is the corresponding additive white Gaussian noise (AWGN). We collect 
all channel vectors into a channel state matrix H(t) = [hi(t), ...,h/^(t)] G C^^^^. Without loss 
of fundamental generality, we assume that the channel coefficients have mean and variance 1 
(e.g., in the case of Rayleigh fading). At the beginning of each slot t, the BS has knowledge 
of the CSI H(t) = [hi(t), h/^ (t)] G C^''^^^^, obtained by some form of training, channel 
prediction and feedback, as discussed in Section HI We assume that H(t) and H(t) are jointly 
stationary and ergodic matrix-valued processes. For convenience, we also assume that H(t) is a 
sufficient statistic for the causal estimation of H(t) from the CSI process {H(t)}. 

While the capacity region of the MIMO-BC in the perfect CSI case (i.e., for H(t) = H(t)) is 
well-known [20], the case of imperfect CSI is still open although outer bounds and achievability 

'The generalization to MIMO-OFDM and frequency selective fading is immediate. 



4 



lower bounds exist. In this work we focus on a simple physical layer signaling scheme based on 
linear precoding and independently generated Gaussian user codes. Nevertheless, the general 
scheduling framework developed in this work can be easily extended to other MU-MIMO 
downlink schemes, such as Tomlinson-Harashima precoding [3], Vector Precoding [21] and 
Dirty-Paper Coding [20], [22]. 

With linear precoding, each A;-th user codeword is a M x T space-time array denoted by 
Ufc(t) = {ufcj(t) : i = 1,...,T}. The signal vector transmitted at symbol interval i of slot 
t is given by Xj(t) = '^he following, we let (H, H) denote a pair of jointly 

distributed random matrices with the same joint distribution of (H(t),H(t)) (independent of 
t by stationarity). A linear precoding signaling scheme is defined as a possibly randomized 
function 7 such that 

7(H) = (5]i(H),...,E^(H),r(H)) 

where 5]fc(H) is the spatial-domain transmit covariance matrix of user k, and r(H) is a transmit 
rate allocation vector. Then, upon observation of the CSI H(t), the signaling scheme 7 chooses 
for each user k a Gaussian generated codebook of rate rfc(H(t)), where the codewords Ufc(t) 
have i.i.d. columns generated according to the Gaussian distribution C!N(0, Sfc(H(t))). We say 
that a scheme 7 is feasible with respect to the power constraint P if XlfcLitr (^Sfc(H)j < P 
with probability 1. The set of all feasible schemes is denoted by r(P). For a given 7 and CSI 
value H, the set of users k such that tr(Sfc(H)) > is called the active set, and will be denoted 
by U^(H). 

The above definition of 7 encompasses in full generality all linear precoding strategies based 
on Gaussian random coding, ranging from beamforming to space-time coding. For later use, we 
recall here two well-known choices for the transmit covariance matrices that will be essential in 
the practical scheduling policy of Section |lVl 

1) A popular choice for MU-MIMO linear precoding consists of computing ZFBF "steering 
vectors" by treating the CSI H as if it was the true channel matrix (see for example [8], [10], 
[23] and references therein). In our notation, this corresponds to choosing an active set 'U-y(H) of 
size not larger than rank(H) and, rank-1 transmit covariance matrices Sfe(H) = Pk^^y^ where 
Pfc > is the transmit power allocated to user k, and Vfc is a unit-length vector obtained by 
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calculating the Moore-Penrose pseudo-inverse 

^^■(IX-^) = H(IX^) (^H"(U^)H(U^))~^ (2) 

of the matrix li{Uy) with columns {hj : j G IX^}, and taking the normalized column of H^(lX-y) 
corresponding to user k. In particular, is orthogonal to all {hj : j G 'U-y(H), j 7^ k}. 
2) At the other extreme of the range of possible linear precoding signaling schemes we find the 
classical space-time coding to a single user [13], [24], [25]. In our notation, this corresponds to 
choosing an active set IX^(H) of size 1 and the transmit covariance matrix Sfe(H) = (P/M)I 
for the only k G 1X^(11). Interestingly, ZFBF serves simultaneously up to M active users, each 
with a rank-1 transmit covariances matrix, while space-time coding serves just one active user, 
with a rank-M transmit covariance matrix. 

For a fixed set of transmit covariance matrices (Si, ... , 'Sr), a linear precoding scheme yields 
a Signal-to-Interference plus Noise Ratio (SINR) at receiver k given by 

SINR.(H, El, . . . , S,) A (3) 

We let Rk{t) denote the effective rate of user k on slot t. In general Rk{t) is different from 
the allocated rate rfc(H(t)) since CSI is not perfect. As far as rate allocation is concerned, we 
consider the following two cases: 

1 ) Outage rates: following standard information theoretic arguments (see [26] and references 
therein), under slot-by-slot coding and decoding, receiver k can reliably decode a rate provided 
that no information-outage occurs, i.e., provided that is smaller than the mutual information 
/fc(H, El, . . . , Y^k) = log (1 + SINRfc (H, El, ... , S^)). As a consequence, for a given signaling 
scheme 7 we define the outage rate as the random variable: 

i?fc(H, 7(H)) = rfc(H) X 1 {rfc(H) < 4(H, Si(H), . . . , I],^(H))} (4) 

where 1{^} is the indicator function of an event A. 

2) Optimistic rates: in this case, we assume that some genie-aided rate adaptation scheme is 
able to achieve an effective instantaneous rate equal to the mutual information: 

RkiH, 7(H)) = 4(H, Si(H), . . . , Sx(H)) (5) 

The system model underlying the outage rate assumption consists of standard ARQ protocol 
that removes Rk{t) = rfc(H(t)) bits/channel use from the transmission buffer of active user k 
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if no outage occurs. The system model underlying the optimistic rate assumption corresponds 
to an idealized fast rate adaptation scheme (see for example [27], [28]). Any practical rate 
adaptation scheme yields performance in between the outage and the optimistic rates defined 
above. Therefore, these two extreme cases are relevant in the sense that they provide upper and 
lower bounds to practical adaptive rate schemes. Under either one of the above assumptions, we 
let the effective rate be Rk{t) = i?fc(H(t), 7(H(t))), given by © or by ©. 

in. Optimal downlink scheduling 

The achievable ergodic rate region % for a given set of feasible physical layer signaling 
schemes, is defined as the closure of the convex hull of all achievable ergodic rate points. Under 
a fixed signaling scheme 7 G r(P), user k is served with a long-term average rate Rk = 
\imt^oo]-J2'T^o^k(ti{T),'y(H.{T))) = E[_Rfc(H, 7(H))], where convergence is with probability 
1 because of ergodicity. Since time-sharing between any set of feasible signaling schemes is also 
a feasible scheme, we have: 

, VA;} (6) 

7er(P) 

where "coh" denotes "closure of the convex hull" and where the expectation is with respect to 
the joint probability distribution of (H, H) and 7 (for randomized signaling schemes). 

We consider an "infinite backlog" situation where all the data to be transmitted are available 
at the BS. The goal of the downlink scheduler is to maximize some concave entry wise non- 
decreasing utility function g{-) of the user individual ergodic rates R = (_Ri, . . . ,-Ri^)E The 
problem that we wish to solve is: 

maximize 5'(R), subject to R G 3^ (7) 

Suppose that the solution R* of (|7]) is found. Then, by definition, there exists a (possibly 
randomized) signaling strategy that achieves R*. A feasible scheduling policy is an algorithm that 
chooses at each time t some physical layer signaling scheme 7 G r(P), based on the history of all 
past transmissions and arrivals, on the observation of the CSI and on the knowledge of the joint 
statistics of all variables in the system. We are interested in finding an explicit scheduling policy 

^By entrywise non-decreasing we mean that for all r G and d G M.^ , g{r) > g{r + 5). Also, recall that a concave 
function is continuous in the interior of its domain. Without loss of generality, we consider g{-) with domain R^. 



3^ = coh IJ |r G K+ : i?fc < E i?fc(H, 7(H)) 
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(denoted for brevity by 7*) that achieves R*. Despite the fact that (|7]) is a convex optimization 
problem, a direct solution is generally overly complicated since !Jl does not admit in general a 
simple characterization. For example, Ji is generally not a polytope, and may be described by 
an uncountable number of linear constraints (supporting hyperplanes). 

Fortunately, we can use the framework of [16] and obtain a dynamic scheduling policy that 
operates arbitrarily closely to the optimal point R*. This is obtained in two steps: first, a dynamic 
scheduling policy that achieves the stability of transmission queues whenever the arrival rates 
are inside is obtained. Then, we build "virtual queues" driven by appropriate "virtual arrival 
processes," such that their arrival rates are as close as desired to the desired rate point R*. Our 
analysis extends the results in [16] to this new context and also provides a new and tighter 
bounding analysis for the queues, particularly for general (possibly negative) concave utilities 
that include the proportional fairness utility. We note that it may be possible to pursue utility 
optimization using the alternative stochastic approximation and fluid transformation approaches 
in [29], [30], [31], [32], although these approaches may not yield explicit queue bounds. Further, 
the stochastic approximation techniques in [29], [30] use an infinite running time average of 
transmission rates, whereas our approach does not require an infinite running time average and 
can thus adapt to system changes. 

A. System stability 

Suppose that the data to be transmitted to user 1,. . . ,K arrive to the BS according to a 
stationary and ergodic vector-valued process A{t) = {Ai{t), . . . ,AK{t)), with rate vector A = 
E[A(t)] (expressed in bit/channel use) and such that < v4fc(t) < Amaxj V t, for some constant 
^max < CO- The BS maintains a transmission queue for each user, and we let Qkit) denote the size 
of the A;-th queue buffer at the beginning of slot t. As described in Section HIl Rk{t) bit/channel 
use are removed from queue k during slot t, i.e., Rk{t) represents the instantaneous "service 
rate" of the A;-th queue. Defining Q(t) = (Qi(t), . . . , (5i^(t)) and R(f) = (i?i(t), . . . , i?i^(t)), 
the queues evolution is described by the stochastic difference equation^ 

Q(t + 1) = max{0, Q{t) - R(t)} + A{t) (8) 

We have the following definition [16]: 



^The function max{-, •} is applied componentwise to vectors. 
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Definition 1: A discrete-time queue Qkit) is strongly stable if lim sup^^^^^ j J2t=o IE[Qa,(t)] < 
oo. The system is strongly stable if all queues A; = 1, . . . , K are strongly stable. 

For convenience, throughout this paper we use the term "stability" to refer to strong sta- 
bility. It can be shown [16] that if Qk{t) is strongly stable and Ak{t) is uniformly bounded 
by a finite constant A^aax, as in our case, then limj^oo Qfc(^)/^ = with probability 1 and 
\imt-^oo'^[Qki't)]/t = 0. These properties are referred to as rate stability and mean-rate stability, 
respectively. In particular, rate stability implies that the time average rate of bits going into the 
queue is equal to the time average rate of bits going out of the queue. 

The system stability region is the the closure of the convex hull of all arrival rate points A for 
which there exists a feasible scheduling policy that achieves system stability [16]. The following 
result yields both the system stability region and the dynamic scheduling policy that stabilizes 
the system for any arrival rate point inside the region: 

Theorem 1: Suppose the arrival vector A(t) is i.i.d. over slots with each entry uniformly 
bounded by some finite constant Amax, and that the joint channel state and CSI pair {H(t), H(t)} 
is i.i.d. over slotsoFor the system defined in SectionUIl the system stability region coincides with 
the ergodic rate region Ji given in Furthermore, for any arrival rate point A in the interior of 
% the system is stabilized by the dynamic scheduling policy 7* defined as follows. For Q G 
and H G C^'^^^, consider the signaling scheme with covariance matrices (S^(H), . . . , S^'(H)) 
and rate allocation vector r*(H) given by the solution of the weighted sum-rate maximization: 



maximize J2k=i Qk^ 



(9) 



-Rfc(H, El, . . . , Ex, r)| H 
subject to ^"^^ tr(Sfe) < P, > V A;, r > 
Then, the dynamic scheduling policy 7* chooses at each time t the signaling scheme defined by 
^ for the current queue states (i.e., for Q = Q(t)) and the current CSI (i.e., for H = H(t)). 
Proof: See Appendix U ■ 
Interestingly, the weighted sum-rate maximization in Q defining 7* involves the conditional 
expected service rates for given CSI: in the absence of perfect CSI the BS schedules the users 
on the basis of the MMSE estimation (conditional mean) of their instantaneous service rates. We 
conclude this section with a note on the optimal rate allocation in the stability policy 7*. Under 



■^This result and the result of Theorem |2] are stated for the i.i.d. case, but they can be extended to jointly ergodic processes 
subject to some mild technical conditions by following the technique of [18]. We omit this extension for brevity. 
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the optimistic rate assumption, r is irrelevant since ([5]) does not depend on r. Under the outage 
rate assumption, using (H]) we obtain the optimal rate allocation r* for a given set of covariance 
matrices and CSI as the solution of (see also [13]): 



rl(H) = argmax r [l - P flog (1 + SINRfe(H, Si, ... , S^)) < r| H 

r>0 L V 



(10) 



B. System optimization 

Going back to the original problem ([7]) and following [16], we build "virtual queues" with 
arrival rate Ji-tuple A arbitrarily close to the desired optimal point R*, although the latter 
is not known a priori. Then, using the stability policy 7* applied to the virtual queues, the 
system necessarily operates at a throughput point R > A (componentwise domination). From 
the monotonicity of g{-) we are guaranteed that the system will operate arbitrarily close to the 
optimal point. 

Specifically, we define 7* as follows: let V, Amax > be suitable constants. At each time t, 
let Akit) = ttk, where a = (ai, . . . , uk) is the solution of 

K 

max Vg{a.) -S^ akQkit) (11) 

a:0<afe<Aniax, Vfc ^-^ 

k=l 

Then, for given Q(t) and CSI H(t) the signaling scheme 7* given in ^ is applied, resulting 
in the service rates R(t). Finally, the virtual queues are updated according to ([8]), with arrivals 
A(t) given by (fTT|) and service rates R(t). The performance of the scheduling policy 7* is given 
by the following: 

Theorem 2: Suppose the joint channel state and CSI pair {H(t),H(t)} is i.i.d over slots. 
Consider the scheduling policy 7* defined above, for given constants V > Q and Amax > 0. 
Assume that g{-) is concave and entry-wise non-decreasing and that there exists at least one 
point r G 31 with strictly positive entries such that g{r/2) > —00. Then: 

(a) The utility associated with the time average transmission rates achieved by 7* satisfies: 

\unMg I I VE[R(r)] I > g{R*{A^^,)) - C/V (12) 

where 

i^(ALx + ]E[log^(l + |h.(t)pP/iVo)]) ^^^^ 
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and where R*(y4jnax) denotes the solution of the problem ^ with the additional constraint 
< ^fc < Amax for all A; = 1, ... , K. 

(b) For any point R G 3i such that < ^a.- < ^max for all k, and for any value /? E [0, 1] we 
have: 

r lV^V-5^r^.M^ g' + V^[^(R^(^max)) - gjm] ,,,, 

hmsup - yyRkE[Qk{T)\ < (14) 

t^oo t f— ' J- — P 

r=U fe=l 

Thus, all queues Qkit) are strongly stable. 

Proof: See Appendix U ■ 

Theorem [21 implies that if Amax is sufficiently large, such that Amax > Rl fo^ ^' then: 

t~i 



liminf 1 VE[R(r)] ) > ^^(R*) - C/V. 



\ T=0 

Hence, the control parameter V can be chosen as large as desired to make the achieved utility 
arbitrarily close to the optimal utility g{R.*) for the problem (|7]). This comes with a tradeoff in 
the virtual queue average sizes that, as seen from (O, grow linearly with V . The virtual queue 
sizes represent the difference between the virtual bits admitted into the queues and the actual 
bits transmitted, and thus affect the time-scales required for the time averages to become close 
to their limiting values. 

C. Proportional fairness and hard fairness scheduling 

We shall focus on two particularly relevant special cases for the system utility function g{-) 
that reflect useful forms of fairness. The proportional fairness schedulling (PFS) is defined by 
the utility function ^(R) = XlfcloS-^fc [5]^ [6]- In this case, the solution of (fTT)) is given by: 

Ak{t) = min j — A^ax 1 (15) 



Qkity 

The hard fairness scheduling (HPS), uses the utility function 5'(R) = min{_Ri, . . . , Rr}- Making 
use of an auxiliary variable p, (fTTI) can be re-stated as 

K 

max - VI akQk{t) (16) 

k=l 

Solving first with respect to a for fixed p and then solving with respect to p we obtain 

iiv>EtMt) ^^^^ 

else 
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An interesting and new aspect of Theorem |2] is that part (b) allows us to find explicitly the tighter 
upper bound on the virtual queue sizes. For example, consider the proportional fairness utility, 
choose a sufficiently large A^ax such that Amax > Rt ^' consider the vector R = R* 

in (fT4l) . Then, we have g{P'W') = K \og{l3)+g{YC) and the bound in Theorem [2] part (b) becomes 
\ where [3* is the unique solution in [0, 1] of the equation log/3 + ^ = 1 + 
For the hard fairness utility, choosing again Amax > R*k for ^ letting R = R* in 
the bound, we have g{[3JV') = [3g*, where g* is the max-min ergodic per-user rate. Hence, the 
bound in Theorem [2] part (b) becomes simply C + Vg*. 

IV. Algorithms 

In order to implement 7*, a simple and easily computable expression for the conditional 
average service rates E[i?fc(H, Si, . . . , S/^, r)|H] is needed. Unfortunately, this is generally not 
available even in the case where the joint statistics of H and H is known and easy to characterize 
(e.g., when H and H are jointly Gaussian). Another difficulty in the implementation of 7* is 
that the weighted sum rate maximization in ^ is generally a non-convex problem for the linear 
precoding signaling schemes at hand. 

To overcome these difficulties, we introduce some approximations. We start by considering 
the ZFBF scheme reviewed in Section HIl In this case, the maximization in ^ reduces to the 
maximization over all active sets U C {1, . . . , K} such that |1X| < M of the weighted sum-rate 

5(lX,p,Q,H) = ^gfcE 

k&U 

where p = {pk} and v^. are the ZFBF vectors defined in Section [III Notice that {v^ : A; G 11} 
and p are functions of H, even though we omit the explicit dependence for the sake of notation 
simplicity. Also, it is understood that = for all k ^U. 

Motivated by the findings on channel state prediction error mentioned in Section H] and 
illustrated in Section |Vl we make the working assumption that the users can be partitioned into 
two classes: the subset IKpr of users with very small channel state prediction error ("predictable 
users") and the subset X^pv of users with very large channel state prediction error ("non- 
predictable users"). In order to develop some intuition, we assume the CSI model = + 
where the CSI estimation error is statistically independent of h^, with mean zero and variance 
CTg per component. Furthermore, we restrict to the on-off power allocation pk = P/\U\ x l{k E 



log 1 + 









No + Y^j(zU:jjtk 







H 



(18) 
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11} , that is known to yield near-optimal sum-rate for the optimal choice of the active user set 
U and sufficiently large P/Nq. We wish to understand whether a given user should be included 
in the active set in the maximization of (fTSi) . For this purpose, we evaluate the gap between the 
actual service rate of user k and the service rate user k would achieve if the BS has perfect 
knowledge of h^. We will evaluate this gap under both assumptions k e Xpj. (corresponding 
to (Tg ^ 0) and k E 3Cnpr (corresponding to cTg ^ 1), and eventually conclude that if a user 
k E %npr is to be served, then no other user should be served in the same slot. 

Suppose that a genie provides the true channel vector to the BS. Then perfect zero-forcing 
to user k is possible. We will denote by H|™''' the genie-aided CSI obtained by replacing 
with hfc in the CSI matrix H. The beamforming vectors computed using the genie-aided CSI 
are denoted by {v|'^'^"^ : j E U}, and have the property that h^v|'^'"'^ = for j 7^ A;, since is 
known perfectly. In general, v|'^'^"^ 7^ unless 11X1 = M. 

The conditional expected service rate under the augmented CSI for user k is given by 



E 



log 1 







2p 




1 H gcnic 


2p 



H,hz 



log 1 



No\U\ 

while in the case of the actual CSI we have 



E 



log 1 + 







2p 






2p 



H 



(19) 



(20) 



Defining the conditional rate-gap as ARk(H.,hk) = -Rfc(H 



gcnic \ 



i?fc(H), using dH and 



the monotonicity of log(-) and Jensen's inequality, after simple algebra we obtain the rate-gap 
upper bound: 

K^k\'P\ , , A , crlm - 1) 



ARkiUM < log 1 



No\U\ 



E 



log 1 + 



No\U\ 



H 



log 1 



iVnlUI 



By the properties of the Moore-Penrose pseudo-inverse Q, we have that E[6fc] > where 
equality holds exactly when |1X| = M and approximately when ^ 0. It follows that if 
k E %pr, then Ai?fc(H, h^) ^ with high probability, i.e., for predictable users the gap between 
perfect and non-perfect CSI is very small, as we may expect. In contrast, if A; G %npr, then 
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Gfc, independent of h^. In this case, Ai?fc(H, h^) ^ + log (l + ^jv'Ll^^ ) grows on 



No\U\ 

average like log(P/Ai'o), unless we let IX = {A;}, i.e., we schedule only user k. 

The above argument leads to the following conclusions: a near-optimal scheduling policy 
should select at each slot t either a group of predictable users and serve them using ZFBF 
spatial multiplexing mode, or a single non-predictable user and serve it using space-time coding, 
that does not require CSI at the transmitter apart from the rate allocation. Operating along these 
guidelines, in all cases the rate-gap with respect to perfect CSI is a constant that does not grow 
with P/Nq. 

For U C 3Cpr, the objective function in (fTSi) becomes 



5P^(1X, p, Q, H) ^ V Q, log ( 1 + ) (21) 



where, for each such subset, the ZFBF vectors {v^} are obtained as the normalized columns of 
the Moore-Penrose pseudo-inverse Then, the power allocation vector p is obtained from the 
standard waterfilling formula [2], [4]. If the number of predictable users is large, the near-optimal 
user selection algorithm of [4] can be used to avoid the combinatorial search over all U C DCpr. 
For the non-predictable users, the corresponding objective function is given by 

-'- + AiNo ) optimistic rates 

Qk maXr>o |r 1 - P (^jh^p < pf^^jj^^ } for outage rates 
Notice that the rate allocation in the outage rate case is actually very simple: it is sufficient to 
know the cdf of jh^p, which is either known a priori or it can be learned "on-line" from the 
channel measurements at each UT. 

The the proposed simplified scheduling policy can be summarized as follows: for a desired 
concave non-decreasing utility function g{-) of the ergodic rates, the virtual arrival processes 
and the corresponding virtual queues are defined in Section Illl-B[ yielding queue buffers Q(t) 
at each scheduling slot t. The scheduler computes S^^^{t) = max^cXpr 'S'^'^'CU, p, Q(t), H(t)) 
and S^^^^{t) = maxfcgXnpr S^'^^{{k}, Q(t)) and chooses to serve the best subset of predictable 
users if SZ.{t) < ^Jnax(^). or the best non-predictable user if S^^it) > S^Ut). 



V. Results and discussion 

In this section we illustrate the performance advantages of the proposed MU-MIMO scheduling 
policies over a conventional "mismatched" PFS scheme that treats the existing CSI as if it was 
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perfect. The mismatched scheme computes the ZFBF vectors {v^} from the CSI matrix H(t) as 
described in Section [III and selects the active user subset by maximizing the mismatched weighted 
sum rate J2k=i''^i'i'^) f^^it), where rk{t) = log ^1 + IMlLzfcLPfe j and where the powers {pk} 
are computed by waterfilling [4]. The coefficients 7k{t) represent time-averaged rates, that are 
updated according to the rule [5], [6] 

7k{t + 1) = (1 - 1/Q7,{t) + {l/QR,{t) 

where Rk{t) denotes the actual service rate of user k, under the outage or optimistic rate 
assumption as defined in & and in ([5]). It is well-known that this algorithm approximately 
maximizes ^i^^ogRk in the case of perfect CSI (i.e., when H(t) = H(t) for all t), when is 
very large. 

In Section IV-Ai we consider an idealized setting where all channels are i.i.d. Rayleigh fading, 
and where a subset !Xpr of users have perfect CSI while the complement set Xnpv has channels 
completely unknown to the BS. Then, in Section IV-BI we considered the SCM channel model 
used in 3GPP standardization [14], and actual channel estimation and prediction schemes. 
Interestingly, even in this very realistic setting, similar performance trends are observed. 

A. Rayleigh fading 

We consider a BS with M = 4 antennas and K = 8 users, with Xnpr = {1)2}. In this section 
we consider the extreme case where the channel of users 1 and 2 is completely unknown to the 
BS, while the channels of the other users are perfectly known. All channel vectors are i.i.d. across 
scheduling slots and in the antenna domain, with elements ~ C!N(0, 1) (independent block-fading 
with spatially white Rayleigh fading). In this case, the mismatched scheme computes the ZFBF 
beamforming vectors for users 3, . . . , 8 without any orthogonality constraint with respect to the 
channels of users 1 and 2, since the latter are unknown and isotropically distributed. 

We start by illustrating the effect of the constants V and Amax on the scheduling performance. 
Fig. \T\ shows the time evolution of the long-term time-average rates achieved by the proposed 
approximation of the policy 7* in the HFS case (max-min throughput), for Amax = 100 and 
V = 100, V = 1000 when SNR is 20dB. In agreement with Theorem |2l by increasing V the time 
response of the algorithm becomes slower while the achieved utility function value improves. In 
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general, the two parameters V and A^ax should be tuned using the bounds of Theorem [2] and 
are functions of channel statistics, of K and M and of the SNR P/Nq. 

Next, we examine the ergodic sum rate and sum log-rate achieved by the new algorithms under 
PFS and HFS, and compare their performance with that of the mismatched PFS scheme. Figs. 
[2] and [3] show the scheduling algorithms performance versus SNR in dB, for both the optimistic 
and the outage rate assumption. The gain of the novel algorithms over mismatched PFS is very 
large, especially under the outage rate assumption. This fact is understood by considering Fig. 
m showing the users "activity fractions", i.e., the fraction of time slots in which a given user 
is active. The mismatched PFS allocates a very large fraction of slots to the non-predictable 
users. This is because if some users have poor quality CSI and the scheduler does not take this 
explicitly into account, then the fairness induced by the PFS utility function forces the system 
to serve these users very often. Hence, the unpredictable users "drain" a large fraction of the 
system capacity despite the fact that there might be a large number of users with very good 
quality CSI. In contrast, the novel schemes treat the non-predictable users separately, and this 
has a very significant impact not only on the ergodic rates of these users, but also on the whole 
system sum rate. It is also interesting to notice that, under the proposed scheduling policies, the 
gap between optimistic rates and outage rates is very small. This indicates that any suitable fast 
rate adaptation (e.g., based on rateless coding and/or incremental redundancy ARQ) has only 
a minor impact on the system performance with respect to a much simpler conventional ARQ 
scheme. 

B. 3GPP channel model and actual channel prediction schemes 

We run extensive experiments based on the so-called "Spatial Channel Model" (SCM) [14]. 
This channel model is not block-fading and the channel coefficients vary continuously over 
time. Although this model is frequency selective, we considered a frequency-flat version of the 
channel corresponding to a single subcarrier of an OFDM system, for consistency with the rest 
of the paper. For a generic user and antenna (indices are omitted), this channel model yields the 
time-varying channel coefficients in the form 



■n 




(23) 



r=l 
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where r] is the number of impinging scattered wavefronts arriving at the receiver (rj = 20 is 
specified in [14]), Aj. are random complex amplitude coefficients, Cr is the Doppler frequency 
shift corresponding to the r-th wavefront, normalized by the signal bandwidth and i ticks at the 
symbol rate. In turns, the Doppler shifts are given by Cr = cos{9r — 6^), where fc is the 
carrier frequency, v is the mobile speed, c denotes light speed, Tg is the symbol interval, 9r is 
the angle of arrival (AoA) of the r-th wavefront, and 9^ is the mobile azimuth direction. 

We assume that a set of M orthogonal downlink pilot symbols are sent by the BS every slot 
of T symbols. Each user estimates and predicts the channel on the next slot using the pilot 
symbols. After thorough comparisons of various channel estimation and prediction schemes, 
not reported here for the sake of space limitation, we report here the results for the two most 
promising schemes in terms of performance versus complexity. The first scheme consists of 
a block-by-block prediction based on the parametric estimation of the parameters {77, Ar, Cr} 
in (I23]) using ESPRIT applied to blocks of ^ 1 pilot symbols, as described in [33]. The 
second scheme is a classical Recursive Least-Squares (RLS), approximating a Wiener MMSE 
predictor for the channel vector sampled at the pilot-insertion rate 1/T [34], [35], [36]. In our 
simulations we considered a system with parameters given in Table H that corresponds to a single 
subcarrier of an OFDM system with 256 subcarriers and bandwidth 256 x 15KHz = 3.84MHz. 
We compared the two prediction methods by considering the four possible different scenarios 
of: 1) High speed (v = 75km/h) vs. low speed {v = 5km/h) mobiles, and 2) well- separated 
and packed AoAs of the impinging wavefronts. A known limitation of any estimator of a linear 
combination of sinusoids in noise (see [37]) is that the estimation error increases sharply when 
the separation between some of the frequency components falls below some minimum resolution 
that depends on the number of pilots N. On the other hand, the RLS prediction error degrades 
as maxr |Cr| is non-negligible with respect to the pilot insertion rate 1/T. It follows that there 
exists a class of channels with both high mobility and clustered AoAs for which all prediction 
methods essentially fail. This corresponds to the "non-predictable" users said before. 

We considered a BS with M = 4 antennas and K = 8 UTs. We report only the results for 
one scenario because of space limitation, but the same trend is observed in a variety of cases 
(see [15]). We consider high-mobility users with ESPRIT parameter estimation/prediction, where 
users 1 and 2 have clustered AoAs (given in Table Hill) and users 3, . . . , 8 have well-separated 
AoAs (given in Table III). The simulation results are obtained by keeping the AoAs fixed, and by 
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averaging with respect to the amplitudes of the SCM model. Fig. [5] shows the average sum-rate 
for the various scheduling algorithms in this case. We notice that the results for these realistic 
channel models and actual channel estimation and prediction schemes are in agreement with 
those for the i.i.d. Rayleigh fading case. 

Appendix I 
Proofs 

Proof of Theorem [B First, we show that any arrival rate for which the system is strongly 
stable must be in Jl. Suppose that for some uniformly bounded i.i.d. process A(t) with rate A, 
there exists a policy that stabilizes the system. Using the queue buffer evolution equation ([8]), 

assuming Q(0) = for simplicity, and summing with respect to r = 0, . . . , t — 1 we obtain: 

t-i t-i 

Q(t)>5^A(r)-5^R(r) (24) 

r=0 T=0 

where R(t) denotes the service rate achieved by the policy. Dividing by t, taking expectations 
and rearranging terms we arrive at: 

l£E[A(r)]<lgE[R(r)] + ^ (25) 

T = T = 

Using the fact that E[A(r)] = A for all r, we see that the left hand side of the above bound is 
equal to A. Since strong stability with a finite Amax implies mean-rate stability [16], it follows 
that E[Q(t)/t] 0, and so the final term in (l25l) converges to the zero vector. Finally, we have 
E[R(r)] G for all r, and hence i T.l=o^[^{^)] ^ ^ for all t (as this is a convex combination 
of the vectors E[R(r)], and 3^ is a convex set). It follows that A is arbitrarily close to a point in 
[R. Because is closed, we conclude that A G 3^. 

Then, in order to show that 7* stabilizes the system for any A in the interior of % we will 
use the Lyapunov drift approach. Let XL(Q) = | J2k=i Ql denote a Lyapunov function defined 
on . The corresponding one-step Lyapunov drift is given by 

A(Q(t)) = E [£(Q(t + 1)) - i:(Q(t))| Q(t)] (26) 

The following result is standard (see [16] and references therein): 
Fact 1: If there exists constants C > and e > such that 

K 

A(Q(t))<C-e^Qfc(t) (27) 

k=l 



18 



then 



t-l K 



c 

e 



limsup^^^E[gfc(r)] < 

^ r=0 k=l 

and hence each queue Qkit) is strongly stable. (> 
In order to show that (ITTI) holds in our case, we use ([8]) and write 

Qk{t + 1)^ < [Qkit)-Rk{t)f + AUt) + 2Ak{t)max{0,Qk{t)-Rkit)} 

< Qu{tf + Rk{tf + Ak{tf - 2Q,{t) [Rk{t) - A,{t)] (28) 

Summing with respect to k and applying conditional expectation E[-|Q(t)] we arrive at 

^ K K 

A(Q(t)) < -J2^[Rkitf + AkitYlQit)] - J2 Qk{t)E [Rk{t) - Ak{t)\Qit)] (29) 

k=l k=l 

Observing that Rk{t) < log(l + \hk{t)\'^P/No), where the latter is the maximum achievable 
instantaneous rate for user k under perfect CSI as if it was alone in the system, it follows that 

if^^lRkitr + A,{trm)] < I (ALx + IE[log'(l + \hk{t)MNo)]) ^ C < oo (30) 
fc=i 

Next, we shall use the following: 

Lemma 1: Let the service rates {Rk{t)} be obtained by the application of the scheduling 
policy 7*. Then, for any R e D^, we have that 



K 



K 



^Qfe(t)E[i?fc(t)|Q(t)] >Y,Qk{t)Rk 



(31) 



k=l 



k=l 



Proof: Notice that is a convex compact region in . For any fixed non-negative weight 
vector Q, the maximum of the linear function J2k=i Qk'^k of r G 3^ is achieved by some 
7 G r(P). Hence, for any R G 3? and weight vector Q(t), there exists 7 G T{P) such that 

K K 

Y,Qk{t)Rk < X;^'i^WlE[^'^(HW,7(H(t)))] 



k=l 



k=l 
K 



Y^Qkim E[i2,(H(t),7(H(t)))|H(t),7] 



k=l 



< E 



K 



5]i,...,SK,r 



k=l 



max ^gfc(t)E i?fc(H(t),Si,...,Sx,r)|H(t) 



K 



E 



k=l 



E 



i?,(H(t),7*(H(t))) H(t) Q(t) 



(32) 
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Since we assumed that the service rates {Rk{t)} are obtained by applying the policy 7*, then, by 



definition, E 



E 



i?fc(H(t),7*(H(t))) H(t) Q(t) =E[i?fc(t)|Q(t)], and the Lemma is proved. 



Now, let A be in the interior of and let the service rates {Rk{t)} be obtained by 7*. Then, 
there exists a e > such that A + el G 3^. Letting R = A + el in (|3TI) and using Lemma [T] we 
have 

K K K 

Y^QkimW) - Mt)mt)] = I]<5fcW(E[i?fc(t)|Q(t)]-Afc)<e^gfc(t) (33) 

k=l k=l k=l 

Using (l30l) and ([33]) in dlH) we find that the condition (l27]) is satisfied under 7*. 

Proof of Theorem IH It is convenient to define the quantities A(t) = 7 X]t=o ^[^(''")] ^'^^ 
= lE^lEfR-l^)]' where A(t) and R(t) are the virtual arrival process and the service 
rate vector induced by policy 7*. We start with a preliminary fact, the proof of which uses the 
general bound (l25l) and the fact that strong stability and uniformly bounded arrival processes 
implies mean-rate stability (i.e., E[Q(t)]/t 0) [16]. 

Fact 2: Suppose queues Q(t) are strongly stable and there is a finite upper bound A^^^ on 
arrivals every slot. If g{-) is a continuous and entry-wise non-decreasing function, then: 



liminf^(A(t)) < liminf ^(R(t)) 
limsup^(A(t)) < g{R*{A^.^^)) 



From (|29|) . (l30l) and Lemma [H we can write 



K 



A(Q(t)) <c-Yl Qk{t)Rk + QkmiMtmit)] 



(34) 
(35) 


(36) 



k=l 



k=l 



where A(Q(t)) is the Lyapunov drift defined in (l26l) . C is the constant given in (l30l) and 
R = (-Ri, . . . , Rk) is any vector in !R. Following the technique of [16], [17], we subtract a term 
related to the utility function from both sides of (|36l) to yield: 



K 



A{Q{t))-VE[g{AmQ{t)] < C-Y,Qkit)Rk 



k=l 



K 



Y,Qk{t)A,{t)-Vg{A{t)) 



k=l 



Qit) 



Note from the definition of 7* that A{t) is chosen for every t to minimize the right hand side 
over all vectors a that satisfy < 0^ < Amax for all k. Let z be any vector in Ji that satisfies 
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< 2fc < ^max for all k. Thus: 

K K 

A{Q{t))-VE[g{AmQit)] < C -Y,Qk{t)Rk + Y,Qk{t)zk-Vg{-L) 

k=l k=l 

Taking expectations of both sides of the above inequality and using the law of iterated expecta- 
tions yields: 

K 

E[i:(Q(t + l))]-E[£(Q(t))]-VE[(7(A(t))] < C-Y,^[Qk{t)]{Rk-Zk)-Vg{7) 

k=l 

For simplicity, assume that Q(0) = 0. The above inequality holds for all t. Summing the above 
over T E {0, . . . ,t — 1}, dividing by t, rearranging terms, and using non-negativity of £(■) gives: 

^ t-l K 

7 E E nQkir)]{Rk -Zk)<C + Vg(A{t)) - Vgiz) (37) 

T=0 k = l 

where we have used Jensen's inequality in the concave function g(-). The above holds for all t, 
all R G and all z G such that < < Amax for all k. Parts (a) and (b) of Theorem [2] are 
proven by plugging different values into (l37l) . We first prove part (b). 

Proof of part (b). Take any point x E Ji such that < < Amax for all k. Choose R = x 
and z = /3x, for any P E [0, 1]. Then from (l37l) we have: 

r=0 k=l ^ 

At this point, we first prove that the queues are strongly stable and then, using Fact [2l we 
obtain part (b). Notice that g{A{t)) < 5'(Amax)> where A^ax is a vector with all entries equal 
to Amax- Using this bound in (l38l) and taking a limsup yields: 

r ^V^V" , M ^ C + ^^?(An.ax) - VgjP^) 

hm sup - > > XkE[Qk{T)\ < (39) 

^ °° T = k = l ^ 

By assumption, there exists at least one point r G that has all positive entries and such that 
g{r/2) > — oo. Choosing /3 = 1/2 and x = r, it follows that the right-end side of (|39| ) is finite 
and hence all queues are strongly stable. 

Because of strong stability and since the arrival processes are uniformly bounded by A^ax < oo 
by construction, we can apply inequality (l35l) of Fact |2] to the right-end side of (l38l) and obtain 
the result of part (b). 

Proof of part (a). We plug R = z = R*(Amax) into (l37l) and obtain: 

g(A{t)) > (7(R^(Anax)) - C/V 

Taking liminf and using (|34l) in Fact [2] yields the result of part (a). 
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Fig. 1. New HFS, ylmax = 100, V = 100 vs. V = 1000. 
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Fig. 2. Ergodic sum rate, Rayleigh fading. 
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Fig. 3. Sum log ergodic rate, Rayleigh fading. 
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Activity fractions (SNR = 20 dB), Outage Rates Assumption 
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Fig. 4. Activity fractions at SNR — 20 dB, outage rate assumption, Rayleigh fading (Black: Mismatclied PFS, Grey: New 
PFS; White: New HFS). 
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TABLE I 

System parameters for simulation 



Description 


Value 


1/Ts, symbol rate 


15KHz 


/c, Carrier frequency 


2.6GHz 


A'', Number of pilot symbols 


200 


T, Pilot symbol spacing 


20 


dmin, scattering distance 


600m 



TABLE II 

Angles of arrival for well-separated case (in radians), 0^, = 4.4780 radians 





6*2 


03 


04 


05 


06 


07 


08 


09 


010 


4.8328 


5.2210 


5.4479 


5.6090 


5.7340 


5.8360 


5.9223 


5.9970 


6.0629 


6.1219 


Oil 


012 




014 


015 


016 


017 


018 


019 


020 


6.1765 


6.2356 


6.3015 


6.3762 


6.4625 


6.5644 


6.6895 


6.8505 


7.0774 


7.4657 



TABLE III 



Angles of arrival for packed case (in radians), 0^ = 0.6939 radians 





02 




04 


O5 


Oe 


O7 


^8 


O9 


^10 


3.7263 


3.6717 


3.7854 


3.6127 


3.8513 


3.5468 


3.9260 


3.4721 


4.0123 


3.3858 


en 






On 


6*15 


eie 


017 


6I18 


6*19 


020 


4.1142 


3.2838 


4.2393 


3.1588 


4.4003 


2.9977 


4.6272 


2.7708 


5.0155 


2.3826 



29 




Fig. 5. Average sum rate, SCM channel model, ESPRIT prediction, optimistic rates. 



